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Abstract 

This paper studies efficient risk-sharing rules for the concave dom- 
inance order. For a univariate risk, it follows from a comonotone 
dominance principle, due to Landsberger and Meilijson [27], that effi- 
ciency is characterized by a comonotonicity condition. The goal of the 
paper is to generalize the comonotone dominance principle as well as 
the equivalence between efficiency and comonotonicity to the multi- 
dimensional case. The multivariate case is more involved (in particular 
because there is no immediate extension of the notion of comonotonic- 
ity), and it is addressed by using techniques from convex duality and 
optimal transportation. 

Keywords: concave order, stochastic dominance, comonotonicity, effi- 
ciency, multivariate risk-sharing. 

1 Introduction 

Motivation. The aim of this paper is to study Pareto efficient allocations 
of risky consumptions of multiple goods in a contingent exchange economy. 
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In this framework, consumption goods are imperfect substitutes, hence con- 
sumption is measured along several different units, instead of being denom- 
inated in one single monetary value. These units can be for instance ma- 
terial consumption and labor, or future consumptions at various subsequent 
dates, or currency units with limited exchangeability. In this setting, risky 
consumption can no longer be represented as a random variable, but as a 
random vector. 

Agents are assumed to have incomplete preferences associated with the 
concave order: a risk (random vector) X is preferred to a risk Y in the con- 
cave order whenever every risk-averse expected utility decision-maker prefers 
X to Y. Again, these preferences form an incomplete order, hence this as- 
sumption (and its empirical content) may appear as relatively weak. We 
shall see, however, that they in fact lead to strong predictions. 

The motivation of the paper is to characterize efficient allocations for 
the concave order on observable data (for instance, insurance contracts). In 
the case of univariate risk, it is known that efficiency for the concave order 
is equivalent to efficiency for some strictly concave expected utility model, 
which in turn yields a tractable characterization of efficiency: comonotonic- 
ity of the allocations. Allocations are comonotone whenever each agent's 
contingent consumption is a nondecreasing function of the aggregate con- 
sumption. Further, a comonotone dominance principle can also be proven: 
if some initial allocation is not comonotone, there is a comonotone alloca- 
tion such that every agent weakly prefers their contingent consumption in 
the new allocation (at least one preferring strictly). Comonotonicity fully 
characterizes efficiency and it is a testable and tractable property. Moreover, 
as a consequence of the comonotone dominance principle, attention may be 
restricted to the set of comonotone allocations, which is convex and almost 
compact. Hence existence results may be obtained for many risk-sharing 
problems (see for instance [26J in the framework of risk measures, or [9] , [10] 
for classes of law invariant and concave utilities). 

Main results. This paper is devoted to the extension of the comonotone 
dominance result and its application to the characterization of efficient allo- 
cations in the multivariate setting. To this end, a definition of comonotonicity 
for the multivariate case is first needed. Roughly speaking, according to the 
definition of multivariate comonotonicity we adopt, an allocation (Xi, X p ) 
(with each Xi being random vectors) is comonotone if it is efficient for some 
strictly concave expected utility model. By first order conditions, this im- 
plies that there is a random vector Z and convex functions (pi such that 
Xi = Vifi(Z), where V<^j is the gradient map of tfi. This is the definition of 
multivariate comonotonicity used by Ekeland, Galichon and Henry in [20] . 
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The comonotone dominance principle is next extended by solving a varia- 
tional problem. More precisely, given an initial allocation and a collection 
of strictly concave utility functions, we maximize the sum of these utilities 
among allocations that dominate the initial allocation. We prove, and this 
is the hard part of the proof, that the corresponding optimal allocation is 
necessarily comonotone. The precise statement of the multivariate comono- 
tone dominance result is, however more complicated than in the univariate 
case since it requires the use of weak closures and a concept slightly stronger 
than strict convexity. This follows from the fact that the set of multivari- 
ate comonotone allocations is neither convex nor compact (even up to con- 
stants), contrary to the univariate case (counterexamples are given). Finally 
one may wonder whether the equivalence between efficient and comonotone 
allocations is preserved in the multivariate case. The answer is yes, up to 
some (interesting) technicalities. Again the precise statement of the result is 
more complicated than in the univariate case. When applied to the univari- 
ate case, our proof of the comonotone dominance result improves upon all 
existing proofs (see [27], [H] and [SU] )• Indeed, it addresses directly the case 
of many agents, it uses neither the discrete case nor a limiting argument, and 
no hypotheses need be made on the aggregate endowment. 

Literature overview. There is a distinguished tradition in modeling 
preferences by concave dominance. Introduced in economics by Rothschild 
and Stiglitz [33], the concave order has then been used in a wide variety of 
economic contexts. To give a few references, let us mention efficiency pricing 
(Peleg and Yaari [32], Chew and Zilcha [12]), measurement of inequality 
(Atkinson [3]), and finance (Dybvig [IE] . Jouini and Kallal [25]). 

In dimension one, the mutuality principle arose in the early work of Borch 
[6J, Arrow [I], [2] and Wilson [38]; see also LeRoy and Werner [28]. Lands- 
berger and Meilijson [27] proved (for two agents and a discrete setting) that 
any allocation of a given aggregate risk is dominated in the sense of concave 
dominance by a comonotone allocation. This comonotone dominance prin- 
ciple has been extended to the continuous case by limiting arguments (see 
[H] and [30]). It implies the comonotonicity of efficient allocations for the 
concave order. The equivalence between comonotonicity and efficiency was 
only proved recently by Dana [13] for the discrete case and by Dana and 
Meilijson [T3] for the continuous case. This equivalence stimulated a line of 
research on comonotonicity in the insurance and finance literature, see for 
instance Jouini and Napp |23j . [23]. On the empirical side, Townsend [37] 
proposed to test whether the mutuality principle holds in three poor villages 
in southern India while Attanasio and Davis ([!]) worked with US labor data. 
The general findings of these empirical studies is that comonotonicity can be 
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usually strongly rejected. A possible explanation of why efficiency is usually 
not observed in the data is that the aforementioned literature only considers 
risk-sharing in the case of one good (monetary consumption) and does not 
take into account the cross-subsidy effects between several risky goods which 
are only imperfect substitutes. Other papers, such as Brown and Matzkin 
(IB]) have tried to test whether observed market data on prices, aggregate 
endowments and individual incomes satisfy the restrictions that are imposed 
by Walrasian equilibrium. In contrast to this approach, we do not assume 
prices to be available to the researcher. 

The notion of multivariate comonotonicity adopted in this paper coin- 
cides (up to some technical details) with the one originally introduced by 
Ekeland, Galichon and Henry in [20] under the name /i-comonotonicity, in 
the context of risk measures. Galichon and Henry use that concept to gen- 
eralize rank-dependent expected utility in [TJ]. Other proposals for multi- 
variate comonotonicity exist and are reviewed e.g. in [33]; however they do 
not seem to be related to efficient risk-sharing. While the results of j2U] 
are strongly related to maximal correlation functionals and to the quadratic 
optimal transportation problem (and in particular Brenier's seminal paper 
[7]), the present approach will rely on a slightly different optimization prob- 
lem that has some familiarities with the multi-marginals optimal transport 
problem of Gangbo and Swi§ch [22]- 

Organization of the paper. The paper is organized as follows. Section 
[2] recalls some definitions and various characterizations of comonotonicity in 
the univariate case. Section [3J revisits the comonotone dominance principle of 
[27] and characterizes efficient risk sharing in the univariate case. A notion 
of multivariate comonotonicity is introduced in Section HI an analogue of 
the comonotone dominance principle is stated, and efficient sharing-rules 
are characterized as the weak closure of comonotone allocations. Section 
concludes the paper. Proofs are gathered in section 

2 Preliminaries 

Given as primitive is a probability space (Q, J 7 , P). For every (univariate or 
multivariate) random vector X on such space, the law of X is denoted C(X). 
Two random vectors X and Y are called equivalent in distribution (denoted 
X ~ Y), if C(X) = C{Y). 

Definition 2.1. Let X and Y be bounded random vectors with values in M. d , 
then X dominates Y for the concave order, denoted X )? Y, if and only if 
K(ip(X)) < K(ip(Y)) for every convex function (f : lR d — > R. If, in addition, 
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K(tp(X)) < K((p(Y)) for some convex function ip, then X is said to dominate 
Y strictly. 

As the paper makes extensive use of convex analysis (Legendre trans- 
forms, infimal convolutions, convex duality), the concave order is defined 
here in terms of convex loss functions while usually defined with concave 
utilities. Clearly the definition above coincides with the standard one. As 
X )p Y implies that E(X) = E(K), comparing risks for ^= only makes sense 
for random vectors with the same mean. We refer to Rothschild and Stiglitz 
[M] and Follmer and Schied [21] for various characterizations of concave 
dominance in the univariate case and to Miiller and Stoyan [31] for the mul- 
tivariate case. Using a classical result of Cartier, Fell and Meyer (see [H] or 
[36]). one deduces a convenient characterization (see section [6] for a proof) of 
strict dominance as follows: 

Lemma 2.2. Let X and Y be bounded random vectors with values in M. d , 
then the following statements are equivalent: 

1. X strictly dominates Y , 

2. X ^Y andC(X) ^ C(Y), 

3. X >^Y and for every strictly convex function ip, K(ip(X)) < K(ip(Y)). 

Given X £ L°°(f2,IR d ), a random vector of aggregate risk of dimension 
d > 1, the set of admissible allocations or risk-sharing of X among p agents 
is denoted A(X): 

p 

A(X) := {Y = (Y u ...,Y P ) G L°°(Q,R d ) : J2 Yi = X ^ 

t=i 

For simplicity, the dependence of A{X) on the number p of agents does not 
appear explicitly. A concept of dominance for allocations of X is defined 
next. 

Definition 2.3. For d > I, let X = (X t , ...,X p ) and Y := (Y u ...,Y p ) be in 
A(X). Then X is said to dominate Y if Xi )p Y± for every i £ {l,...,p}. 
If, in addition, there is an i £ {1, such that X{ strictly dominates Y i} 

then X is said to strictly dominate Y . An allocation X £ A(X) is Pareto- 
efficient (for the concave order) if there is no allocation in A(X) that strictly 
dominates X. 
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It may easily be verified that dominance of allocations can also be defined 
as follows. Let X and Y be in A(X), then X dominates Y if and only if 

p v 
E(5>,(JQ)<E(^(Y,)), (2.1) 
i=i i=i 

for every collection of convex functions (pi : M. d — > R. Moreover, X strictly 
dominates Y if and only if the previous inequality is strict for some collection 
of convex functions (pi : R d — > R . Note that from lemma 12.2} it is equivalent 
to require that the inequality is strict for every collection of strictly convex 
functions. Therefore, if X is the solution of the problem 

v 

inf{ j^EfopS)) : (Y 1 ,...,Y p )eA(X)^ (2.2) 

i=l 

for some collection of strictly convex functions tp i} then X is efficient. Finally, 
recall that in the univariate case, comonotonicity is defined by: 

Definition 2.4. A collection (Xi, X p ) of p real-valued random variables 
on (flyJ 7 , P) is comonotone if for every G {l,...,p} 2 , 

(Xi(u') -JSQHXXjV) -Xj(u))) > for¥®F-a.e. (u,u/) G fi 2 . 

It is well-known that comonotonicity of (Xi, ...,X P ) is equivalent to the 
fact that each Xi can be written as a nondecreasing function of the sum 
X = (see for instance Denneberg [15]). Therefore (Xi, X p ) is 

comonotone if and only if there are nondecreasing functions fi summing 
to the identity such that Xi = fi(X). Note that the functions fi are all 
1-Lipschitz . The extension of this notion to the multivariate case (i.e when 
each Xi is Revalued) is not immediately obvious and will be addressed in 
Section HI 

We now provide another characterization of comonotonicity based on the 
notion of maximal correlation. From now on, assume that the underlying 
probability space (fi, J 7 , IP) is non-atomic which means that there is no A G J 7 
such that for every B G 7 if ¥(B) < F(A) then ¥(B) = 0. It is well-known 
that (fi, J 7 , P) is non-atomic if and only if a random variable that is uniformly 
distributed on [0, 1], which is denoted U ~ li ([0, 1]), can be constructed on 
(fi, ^P). Let Z G L\n, J^P), and define for every X G L°°(fi, J^P), both 
Z and X being univariate here, the maximal correlation functional: 

q z (X) := sup E(ZX) = supE(ZX) = sup E(ZX). (2.3) 

Xr^X Zr^Z Zr^Z, X ^ X 
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The functional Qz has extensively been discussed in economics and in finance, 
therefore only a few useful facts are recalled. Let F^ 1 be the quantile function 
of X, that is the pseudo- inverse of distribution function Fx- From Hardy- 
Littlewood's inequality, one has 

qz(X)= f F x \t)F z \t)dt, 
Jo 

and the supremum in (12. 3p is achieved by any pair (Z, X) of comonotone ran- 
dom variables (F^ 1 (U) , F x x (U)) for U uniformly distributed. By symmetry, 
one can either fix Z or fix X. Fixing for instance Z, the supremum is achieved 
by F X \U) where U ~ U ([0, 1]) and Z = F z l {U). When Z is non-atomic, 
there is a unique U = F Z (Z) such that Z = F^^U), and the supremum is 
uniquely attained by the non-decreasing function of Z, F x o Fz(Z): 

qz(X) = E(ZF x 1 oF z (Z)). (2.4) 

Also note that Qz is subadditive: Qz(52i Xi) < J2i Qz(Xi). 

Proposition 2.5. Let (Xi, X p ) be in L°°(0, J 7 , P). The following asser- 
tions are equivalent: 



1. (Xi,...,X p ) are comonotone, 

2. for any non-atomic Z £ L l (Q, J 7 , P), 

Qz = X>M> (2.5) 

3. for some non-atomic Z £ L X (Q, J 7 , P), i2.5\) holds true. 

Proof. For the sake of simplicity, we restrict ourselves to p = 2 and set 
(X U X 2 ) = (X,Y). Point 1 implies point 2 since F x l +Y = F x l + F" 1 for 
comonotone X and Y. To show that point 3 implies point 1, assume that 
for some non-atomic Z, one has (12.51) . which by sublinearity is equivalent to 
Qz(X + Y) > Qz{X) + Qz(Y)- Let Zx+y (resp. Zx and Zy) be distributed 
as Z and solve sup^ z EiZX) (resp. £z(X) and Qz(Y)). One then has: 

E(Z x+y (X + Y)) > E{Z X X) + E(Z Y Y). 

AsE(Z x+Y X) < E(Z X X) andE(Z x+Y Y) < E(Z Y Y), it follows E(Z X+Y X) = 
E{Z X X) = q z (X) and E(Z X+Y Y) = E(Z Y Y) = g z (Y), hence from (j2~t) . 
X = i 7 ^ 1 o F Zx+y (Z x+ y)) and F = Fy 1 o F Zx+y (Z x +y)), proving comono- 
tonicity. □ 
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Proposition 12.51 was the starting point of Ekeland, Galichon and Henry 
[2"U] for providing a multivariate generalization of the concept of comono- 
tonicity. In the sequel we shall further discuss this multivariate extension 
and compare it with the one proposed in the present paper. 

3 The univariate case 

A landmark result, due to Landsberger and Meilijson j2jj states that any 
allocation is dominated by a comonotone one. The original proof was given 
in the discrete case for two agents, and the results were extended to the 
general case by approximation. We give an alternative proof in the Appendix 
based on the same approach we shall use in the multidimensional case. This 
proof is based on a certain optimization problem; we believe that, even in 
the unidimensional case, it is of interest per se since it does not require 
approximation arguments and slightly improves on the original statement 
by proving strict dominance of non-comonotone allocations. Contrary to 
Landsberger and Meilijson, one needs however to assume, as before, that the 
probability space J 7 , P) is non-atomic. 

Theorem 3.1. Let X be a bounded real-valued random variable on the non- 
atomic probability space (f^J 7 , P), and let X = (A 1; X p ) £ A(X) be an 
allocation. There exists a comonotone allocation in A(X) that dominates X. 
Moreover, i/X is not comonotone, then there exists an allocation that strictly 
dominates X. 

As an application, we have: 

Theorem 3.2. Let X be a bounded real-valued random variable on the non- 
atomic probability space J 7 , P) and let X = (X\, X p ) £ A(X). Then 
the following statements are equivalent: 

1. X is efficient, 

2. X is comonotone, 

3. there exist continuous and strictly convex functions (ipx, ip p ) such 
that X solves 

p p 
inf^E^y,)) : ]Tr 4 = x}. 

1=1 1=1 

Proof. Point 1 implies point 2: the comonotonicity of efficient allocations 
of X follows directly from Theorem 13.11 Point 2 implies point 3: if X = 
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(X 1; ...,X p ) is comonotone, let us write X^ = fi{X) for some nondecreasing 
and 1-Lipschitz functions ff [m, M] —> R (with M := EsssupX, m := 
EssinfX) summing up to the identity map. Extending the fi functions by 
f t (x) = fi(M) + (x - M)/p for x > M and fi(x) = fi{m) + (x - m)/p 
for x < m, one gets 1-Lipschitz nondecreasing functions summing up to the 
identity everywhere. Let <p(x) := J Q X fi(s)ds for every x. The functions (fi 
are convex and C 1 ' 1 (i.e. C 1 with a Lipschitz continuous derivative) and 
have quadratic growth at oo. The convex conjugates^ ipi ■= are strictly 
convex and continuous functions, and by construction, one has for every i, 
X G dtpi(Xi) a.s., which implies that (X\, X p ) minimizes E(^« V'i(^)) 
subject to J2i Yi = X, which proves point 3. Point 3 implies point 1 since 
the ipi functions are strictly convex; if (Xi, ...,X p ) satisfies point 3 then it is 
an efficient allocation of X. 

□ 

Corollary 3.3. Let (fi, J 7 , P) be non-atomic, then the set of efficient alloca- 
tions of X is convex and compact in L°° up to zero-sum translations (which 
means that it can be written as {(Ai, A p ) : X)£=i Ai = 0} + A with A 
compact in L°°). In particular, the set of efficient allocations of X is closed 
in L°°. 

Proof. Let M := EsssupX, m := EssinfX and define Kq as the set of 
functions (/i,...,/ p ) 6 C([m, M], W) such that for each nondecreasing /j, 
fi{ m ) — rajp and Y7i=\ fi( x ) = x f° r ever y % £ [ m , M], and let 

v 

X:=X + {(A 1 ,...,A P ) : X> = °}- 

i=l 

The convexity claim thus follows from theorem 13.21 and the convexity of K. 
Let us remark that elements of Kq have 1-Lipschitz components and are 
bounded. The compactness of K in C([m, M], M?) then follows from Ascoli's 
theorem. The compactness and closedness claims directly follow. □ 

Convexity and compactness of efficient allocations are quite remarkable 
features and as will be shown later, they are no longer true in the multivariate 
case. Note also that efficient allocations are regular: they are 1-Lipschitz 
functions of aggregate risk. 

^Xet us recall that the Legendre transform or convex conjugate of tpi is by definition 
given by (p*(x) := sup {x ■ y ~ <Pi(y)}. 
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4 The multivariate case 



The aim of this section is to generalize to the multivariate case the results 
obtained in the univariate case. More particularly, Landsberger and Meil- 
ij son's comonotone dominance principle are extended: 1) any allocation is 
dominated by a comonotone allocation; 2) any non comonotone allocation is 
strictly dominated by a comonotone one. 

When addressing these generalizations it is not immediately clear what 
is the appropriate notion of comonotonicity in the multivariate framework. 
Let us informally give an intuitive presentation of the approach developed 
in the following paragraphs. A natural generalization of monotone maps in 
several dimensions is given by subgradients of convex functions. It is there- 
fore tempting to say that an allocation (X\, ...,X P ) e -A(X) is comonotone 
whenever there is a common random vector Z (interpreted as a price) and 
convex functions tfi (interpreted as individual costs) such that Xi e dtfi(Z) 
a.s. for every i. Formally, this is nothing but the optimality condition for 
the risk-sharing or infimal convolution problem 



where ipi = (p* (the Legendre Transform of tpi). This suggests a definition of 
comonotone allocations as the allocations that solve a risk-sharing problem of 
the type above. This has a natural interpretation in terms of risk-sharing, but 
one has to be cautious about such a definition whenever the ipi functions are 
degenerat^]. Indeed, if all the ^ functions are constant, then any allocation 
is comonotone in that sense! This means that one has to impose strict con- 
vexity in the definition. We shall actually go one step further in quantifying 
strict convexity as follows. Given an arbitrary collection w = (wi, ...,w p ) of 
strictly convex functions, we will say that an allocation is w-strictly comono- 
tone whenever it solves a risk-sharing problem of the form (14.1 1) for some ipi 
functions which are more convex than the Wi (i.e. ipi — W{ is convex for every 
i). Allocations which can be approached (in law) by strictly u>-comonotone 
will be called comonotone. Since they solve a strictly convex risk-sharing 
problem, w-strictly comonotone allocations are efficient and the main goal 
of this section will be to generalize the univariate comonotone dominance 
result. We shall indeed prove that for any allocation X e -A(X) and any 

2 In the univariate case, the situation is much simpler since one can take Z = X, and 
since the Xi variables sum up to X, each convex function ipi has to be differcntiablc i.e. 
all the tpi necessarily are strictly convex. In other words, degeneracies can be ruled out 
easily in the univariate case. 



p 
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choice of w, there is a u>-comonotone allocation Y e A(X) that dominates 
X (strictly whenever X is not itself w-comonotone) . The full proof is detailed 
in Section El but its starting point is quite intuitive and consists of studying 
the optimization problem: 

v 

infj^EK^)) : {Y 1 ,...,Y p )eA{X), Y % ^ X u z = l,...,pj. (4.2) 

i=l 

Clearly, the solution Y of (14. 2 p dominates X. A careful study of the dual of 
(14.21) will enable us to prove that Y is necessarily u>-comonotone, thus giving 
the desired multivariate extension of Landsberger and Meilijson's comono- 
tone dominance principle. Note also, that our proof is constructive since it 
relies on an explicit (although difficult to solve in practice) convex minimiza- 
tion problem. 

This section is organized as follows. In paragraph 14.11 we shall refor- 
mulate the problem in terms of joint laws rather than random allocations. 
This is purely technical but will enable us to gain some linearity and some 
compactness in (14.21) . We then define precisely our concepts of multivari- 
ate comonotonicity in paragraph 14.21 Paragraph 14.31 states the multivariate 
comonotone dominance result, i.e. the multivariate generalization of Lands- 
berger and Meilijson's results. Finally, in paragraph 14.41 we gather several 
remarks on multivariate comonotonicity and emphasize some important qual- 
itative differences between the univariate and multivariate cases. 



4.1 From random vectors to joint laws 

From now on, it is assumed that the underlying probability space (Q, J 7 , P) 
is non-atomic, that there are p agents and that risk is <i-dimensional. X is 
a given R^-valued L°° random vector modeling an aggregate random multi- 
variate risk, while X = (X\, X p ) is a given L°° sharing of X among the p 
agents, that is 

i=i 

Let 70 := £(X) be the joint law of X and mo := C{X). Let 7 be a probability 
measure on (§L d ) p and 7* denote its i-th marginal. Note that, C(Yi) is the 
i-th marginal of C(Y). Let IIe7 be the probability measure on M. d defined 
by 

f ip(z)dU sl (z)= [ ^V^,..,^^^^!), (4.3) 
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(where Cq denotes the space of continuous functions that tend to at do). It 
follows from this definition that if 7 = £(Y), then Hs7 = £(X^)- Hence, 
if Y G A(X) and 7 = £(Y), then IIs7 = m = C{X). In other words, if 
7 = £(Y) with Y G A{X), then 

J <p( Xl + ... + x d )d-f(x u x d ) = J ip(z)dm (z), G C (R d , E). (4.4) 

Since Y is bounded, 7 is compactly supported. It follows from the next 
lemma that {£(Y), Y G *4(X)} coincides with the set of compactly sup- 
ported probability measures 7 on (IR d ) p that satisfy (I4.4p : 

Lemma 4.1. Assume (Q, J 7 , P) is non-atomic. If 7 is a compactly supported 
probability measure on (M. d ) p and satisfies l[4-4\ ), then there exists a random 
vector Y = (Y U ...,Y P ) G A(X) such that £(Y) = 7. Hence {£(Y), Y G 
A(X)} = JH(m ), where Ai(m ) is the set of compactly supported probability 
measures on (IR d ) p such that He7 = m = IIs7o- 

In the sequel, joint laws M.(mo) will be used instead of admissible allo- 
cations A(X). For compactness issues, a closed ball B G lR d centered at 
such that m is supported by B p is chosen, and attention is restricted to the 
set of elements of M.(mo) supported by pB (meaning that only risk-sharings 
of X whose components take value in B will be considered). We thus define 

M B (m ) := {7 G M(m ) : 7 (5 P ) = 1}. 



4.2 Efficiency and comonotonicity in the multivariate 

case 

Let C be the cone of convex and continuous functions on B, dominance and 
efficiency in terms of joint laws are defined as follows: 

Definition 4.2. Let 7 and ir be in M.B{ m o), then 7 dominates tt whenever 
C v c v 

\ y^ j ip i (x i )d'y(xi, ...,x p ) < / y2(pi(x i )dir(x 1 ,...,x p ) (4.5) 

J BP i=1 Jbp - =1 

for all functions (<fi, <p p ) G C p . If, in addition, inequality \4-5\) ^ s strict 
whenever the tpi functions are further assumed to be strictly convex, then 7 is 
said to dominate strictly it. The allocation 7 G A^s(m ) is efficient if there 
is no other allocation in MB(mo) that strictly dominates it. 
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Given 7 G Mb (mo), it is easy to check (taking functions <fi(x) = \xi\ n 
in ( 14 .5p and letting n — > oo) that any 7 G .M(m ) dominating 70 (without 
the restriction that it is supported on B p ) actually belongs to .Ms (mo). 
Hence the choice to only consider allocations supported by B p is in fact not 
restrictive. Indeed, if 7 is supported by B p , then efficiency of 7 in the usual 
sense, i.e. without restricting to competitors supported by B p , is equivalent 
to efficiency among competitors supported by B p . 

To define comonotonicity, let ip := (ipi, ...,ip p ) be a family of strictly 
convex continuous functions (defined on B). For any x G pB, let us consider 
the risk sharing (or infimal convolution) problem: 

{p p 
^2^i(yi) : y% G B, ^2y i = x 
1=1 1=1 

This problem admits a unique solution which will be denoted 

T 4> (x) :=(T;(x),...,T^)). 

Note that, by definition 

p 

^T^(x) = x, Va; G pB. (4.6) 
i=i 

The map x 1— > T^(x) gives the optimal way to share x so as to minimize 
the total cost when each individual cost is ipi. It defines the efficient allocation 
Tip(X) := (T^(X), ...,T^(X)) with joint law 7^ defined by: 



f{y x ,...,y p )d^(y) := / f(T^,(x))dm (x) 

Bp JpB 

for any / G C(B P ). One then defines comonotonicity as follows: 

Definition 4.3. An allocation 7 G Mg(m ) strictly comonotone if there 
exists a family ip := (ipi, ...,ip p ) of strictly convex continuous functions such 
that 7 = 7^. Given a family w := (wi, ...,w p ) of strictly convex functions in 
C 1 (fi), an allocation 7 G A^B(m ) «s w-strictly comonotone if there exists a 
family ip := ...,ip P ) of convex continuous functions such that ipi — uii G C 
for every i and 7 = 7^. 

We shall soon show that strictly comonotone random vectors are in the 
image of monotone operators (subgradients of convex functions), evaluated 
at the same random vector, p(X), which justifies the terminology "comono- 
tonicity" in the multivariate setting. By definition, any strictly comonotone 
allocation is efficient. As the set of strictly comonotone allocations is not 
closed, we are led to introduce another definition. 
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Definition 4.4. An allocation 7 £ .Ms (mo) is comonotone if there exists 
a sequence of strictly comonotone allocations that weakly star converges to 
7. Given a family w := (wx, w p ) of strictly convex functions in C l (B), 
an allocation 7 £ .Ms (mo) is w- comonotone, if there exists a sequence of 
w-strictly comonotone allocations that weakly star converges to 7. 

Definitions 14.31 and 14.41 will be discussed in more detail in paragraph 14.41 
To understand the previous notions of comonotonicity and in particular why 
these allocations are called comonotone, it is important to understand the 
structure of the maps. 

Let us first ignore regularity issues and further assume that the ^ func- 
tions are smooth as well as 1/)* their Legendre transforms. Without the 
constraints Xi £ B, then the optimality conditions imply that there is some 
multiplier p = p(x) such that 

V^(T;(x))= P , hence, T*(x) = V#(p). 

Using (14.61) . one gets 

p p 
i=i 3=1 

thus, 

t;(x) = vv? |v£ot j . 

The maps Tl are therefore composed of gradients of convex functions that 
sum up to the identity. In dimension 1, gradients of convex functions are 
simply monotone maps (and so are composed of such maps). In higher 
dimensions, a richer and more complicated structure emerges that will be 
discussed later. Let us now consider the full problem with the constraints 
that Xi £ B and still assume that the ipi functions are smooth, then the 
optimality conditions read as the existence of a p and a Aj > such that 
V^i{T^(x)) = p — XiT^(x) holds together with the complementary slackness 
conditions: Aj = whenever Tl{x) lies in the interior of B. 

4.3 A multivariate dominance result and equivalence 
between efficiency and comonotonicity 

Let us fix an allocation X = (Xi, X p ) £ A(X) such that X £ B p a.s., 
and set 70 = £(X) so that 70 £ A^s(m ). A family w := (wx, ■■■,w p ) of C 1 
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functions is also given, each of them being strictly convex on B as in section 
14.21 The first main result in the multivariate case is a dominance result, it 
states that every allocation is dominated by a w-comonotone one and that 
the dominance is strict if the initial allocation is not itself w-comonotone. 

Theorem 4.5. Let 7 = £(X) and w be as above. Then there exists some 
7 G M.(mo) that is w-comonotone and dominates 70. Moreover if 70 is not 
itself w-comonotone, then 7 strictly dominates 70. 

The proof of this result will be given in section [6j Without giving details 
at this point, let us explain the main arguments of the proof: 

• The optimization problem (14.21) admits a unique solution Y with law 
7 = C(Y), which is efficient and dominates 70 = £(X). 

• One then proves that 7 is necessarily w-comonotone, by showing that 
that w-comonotonicity is an optimality condition for (14.21) . As usual in 
convex programming, optimality conditions can be obtained by duality. 
This leads to consider the problem 

v v 
inf {l( " n Mj2 Xi ">) : *Pi- w i convex, V*}. (4.7) 

i=l i=l 

By a careful study of (14 .7p . one can prove (but this is rather technical) 
that 7 is w-comonotone. 

• It remains to show that 7 strictly dominates 70 unless 70 is itself w- 
comonotone. From lemma 12. 2\ it suffices to show that Y ^ X. But if 
70 is not w-comonotone, then X cannot be optimal for (14. 2 j) and thus 
Y ^ X. 

In terms of efficiency, the following thus holds: 
Theorem 4.6. Let 7 G Ai 3(1^0) and w be as before. Then 

1. if '7 is strictly w-comonotone, then it is efficient, 

2. if 7 is efficient, then it is w-comonotone for any w, 

3. the closure for the weak-star topology of efficient allocations coincides 
with the set of w-comonotone allocations (which is therefore indepen- 
dent ofw). 
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Proof. Point 1 is a property already mentioned several times. Point 2 follows 
from Theorem 14.51 and point 3 follows from points 1 and 2. 

□ 

Note that by definition, if 70 is strictly u>-comonotone then the value 
of problem (14. 7\i is zero. We shall also prove (see section E]) weak form of 
the converse, namely that if the value of problem f)4.7p is zero then 70 is 
w-comonotone. Therefore, the value of (14. 7p as a function of the joint law 
70 can be viewed as a numerical criterion for comonotonicity and thus for 
efficiency One can therefore, in principle, use on data this value as a test 
statistic for efficiency. 

4.4 Remarks on multivariate comonotonicity 

Comparison with the notion of /x-comonotonicity of [20] . The notion 
of multivariate comonotonicity considered in this paper is to be related to 
the notion of /i-comonotonicity proposed by Ekeland, Galichon and Henry in 
|2Uj . Recall the alternative characterization of comonotonicity given in the 
univariate case in Proposition 12. 5t X\ and X 2 are comonotone if and only if 
(Xi + X 2 ) = Qfj, (X\) + (X 2 ) for a measure /z that is sufficiently regular. 
In dimension d, [20] have introduced the concept of /i-comonotonicity, based 
on this idea: if fi is a probability measure on M. d which does not give positive 
mass to small sets, two random vectors X\ and X 2 on M d are called fi- 
comonotone if and only if 

q„ (X 1 + X 2 ) = 6ll (Xx) + £ M (X 2 ) , 

where the (multivariate) maximum correlation functional (see e.g. [35] or 
[20]) is defined by 

6lx (X) = sup e(x-Y\. 

The authors of [2U] show that X± and X 2 are /x-comonotone if and only if 
there are two convex functions ipi and ip 2 , and a random vector U ~ /i such 
that 

X x = Vipx (U) and X 2 = V^ 2 (U) 

holds almost surely. Therefore, the present notion of multivariate comono- 
tonicity approximately consists of calling Xi and X 2 comonotone if and only 
if there is some measure fi such that X\ and X 2 are /i-comonotone. There 
are, however, qualifications to be added. Indeed, [2U] require some regular- 
ity on the measure \i. In the current setting, no regularity restrictions are 
imposed on //; but instead restrictions on the convexity of and %fj 2 have to 
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be imposed to define the notion of w-comonotonicity before passing to the 
limit. Although not equivalent, these two sets of restrictions originate from 
the same concern: two random vectors are always optimally coupled with 
very degenerate distributions, such as the distribution of constant vectors. 
Therefore one needs to exclude these degenerate cases in order to avoid a 
definition which would be void of substance. This is the very reason why the 
strictly convex Wi functions had to be introduced. 

Comonotone allocations do not form a bounded set. In the scalar 
case, comonotone allocations are parameterized by the set of nondecreasing 
functions summing to the identity map. This set of functions is convex and 
equilipschitz hence compact (up to adding constants summing up to zero). 
This compactness is no longer true in higher dimensions (at least when w = 
and we work on the whole space instead of B), and we believe that this is a 
major structural difference with respect to the univariate case. For simplicity 
assume that p = 2. As outlined in paragraph 14.21 a comonotone allocation 
(Xi, X 2 ) of X is given by a pair of functions that are composed of gradients 
of convex functions and sum up to the identity map. It is no longer true, in 
dimension 2 that this set of maps is compact (up to constants). Indeed, let 
us take neN*, and quadratic ipi and ip2 of the form 



Then the corresponding map is linear, and Tl is given by the matrix 



which is unbounded. 

Comonotone allocations do not form a convex set. Another differ- 
ence with the univariate case is that the set of maps of the form used to 
define comonotonicity is not convex. To see this (again in the case p = d = 2), 
it is enough to show that the set of pairs of 2 x 2 matrices 

K := (Si(Si+S2)~ , ^(Si+S^) -1 ), Si symmetric, positive definite, i = 1,2} 

is not convex. First let us remark that if (M 1; M 2 ) G K then M\ and M 2 
have a positive determinant. Now for n G N*, and e G (0, 1) consider 




with 





-i 
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S f _ ( i V n - £ \ s' = ( 1 ~ y/n ~ £ \ 

1 y y/n — e n J ' 2 V — V n ~ £ n J ' 

and define: 

Mi = Si(Si + 5' 2 )~ 1 , M'^S'iW + S'J- 1 , i = l,2. 
If K was convex then the following matrix would have a positive determinant 

/ 1 + YjEI \ 

Mi + M\ = \ ,r— , 2 2n , 

\ 2 "~r 2 1 / 

which is obviously false for n large enough and e small enough. 



5 Concluding remarks 

In this paper, we have extended Landsberger and Meilijson's comonotone 
dominance principle to the multivariate case by introducing the variational 
problem (14.21) . We have then extended the univariate theory of efficient risk- 
sharing to the case of several goods without perfect substitutability, and we 
derived tractable implications. Two observations can be made at this point. 
In the first place, this paper demonstrates the intrinsic difficulty of the multi- 
variate case, as many features of the univariate case do not extend to higher 
dimensions: computational ease, the compactness and convexity of efficient 
risk-sharing allocations. Second, it illustrates the need for qualification in- 
herent to the multivariate case. Contrary to the univariate case, the need to 
quantify strict convexity as in this paper comes by no coincidence. In fact, 
just as [2U] impose regularity conditions on their "baseline measure" to avoid 
degeneracy, we work with cones which are strictly included in the cone of 
convex functions by quantifying the strict convexity of the functions used. 

Getting back to our initial motivation, namely, finding testable implica- 
tions of efficiency for the concave order, we already emphasized in paragraph 
14.31 that one obtains as a byproduct of our variational approach a numerical 
criterion that could in principle be used as a test statistic for comonotonic- 
ity and thus for efficiency. We thus believe that the present work paves the 
way for an interesting research agenda. First of all, an efficient algorithm 
to decide whether a given allocation in the multivariate case is comonotone 
or not remains to be discovered - we are currently investigating this point. 
The convex nature of the underlying optimization problem helps, but the 
constraints of problem (V*) are delicate to handle numerically. Finally, this 
work opens a research agenda on the empirical relevance of the multivari- 
ate theory confronted to the data: do observations of realized allocations of 
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risk satisfy restrictions imposed by multivariate comonotonicity? As men- 
tioned above, tests in the univariate case have been performed by [I] and [37] 
and suggest rejection. But there is hope that in the more flexible setting of 
multivariate risks, efficiency would be less strongly rejected. 

6 Proofs 

6.1 Proof of Lemma 12.21 

Clearly 1 2 and 3 =>- 1 are obvious. To prove that 2 =>■ 3, assume that 2 
holds true. Let \i := C(X) and v := C{Y). These probability measures are 
supported by some closed ball B, and the Cartier- Fell-Meyer theorem states 
that there is a measurable family of conditional probability measures (T x ) x£ b 
such that T x has mean x and for every / continuous function, one has 



Since fi ^ u, fi({x G B : T x ^ 5 X }) > 0, one deduces from Jensen's 
inequality that for every strictly convex function ip, K(if(Y)) > E,(if(X)). 



For notational simplicity, assume that d — 1, p — 2, X takes values in [0,2] 
a.s. (so that m has support in [0, 2]) and 7 is supported by [0, l] 2 . For every 
neW and k e {0, 2 n+1 }, set 



EC/QO) = / f(v)Mv) =11 f(y)dT x ( y )d»(x) 



J B J B J B 



6.2 Proof of Lemma 14.11 



X 



n 




where A 



k,n ■ 




k=0 



and 




Decompose the strip Cfc jn into a partition by triangles 




T, 




i i + 1 r 



\3 j + ir 



'71 ' 2 n 



x 



k<i+j<k+l 



Since lis (7) = ^0 one has: 



P(4k,n) = 7(C fc ,„) = 



k<i+j<k+l 
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and since (O, J 7 , P) is non-atomic, it follows from Lyapunov's convexity the- 
orem (see [22]) that there exists a partition of Ak n into measurable subsets 
A l jf n such that 

7(T^)=P(4y,V(i,j)e{0,...,2"} : fc<i + j<* + l. (6.1) 
Choose (yi,y 2 )fcf n & T^ n and define 

2 n+l 

fc=0 k<i+j<k+l 

We may also choose inductively the partition of Ak, n by the A]f n to be finer 
and finer with respect to n. By construction, one obtains 

max (\\X n - X\\ L ~, \\X n - Y» - Y 2 n \\ L ~, \\Y n+1 - Y n \\ L ^ < 1 

so that Y n is a Cauchy sequence in L°°, and thus converges to some Y = 
(Y"i, Y 2 ). One then sees that Y x + Y 2 = X, and passing to the limit in (16.11) . 
it follows that £(Y) = 7. 

6.3 Proofs and variational characterization for the mul- 
tivariate dominance result 

The proofs will very much rely on the linear programming problem: 
(V*) sup - / ^2w i (x i )d'j(x) 

76-^(70) J BP i=1 

where -^(70) consists of all 7 G .Ms (mo) such that for each i the marginal 
7* dominates the corresponding marginal of 70 i.e.: 

/ (p(xi)d n f(x) < / (p(x i )d / y (x), Vy? convex on B. 
J bp J Bp 

Problem (V*) presents similarities with the problem solved in [22]. In the 
optimal transport problem considered in [22J, one minimizes the average of 
some quadratic function over joint measures having prescribed marginals 
whereas (V*) includes dominance constraints on the marginals. To shorten 
notations, define 

p 

rj(x) := - ^2wi(xi),Wx = (x 1 ,...,x p ) e B p 
1=1 
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(V*) is the dual problem (see the next lemma for details) of 

Up p 
Jp i=i i=\ 

where E consists of all families (p := (ipi, ....,cp p ,if ) G C(B) P x C(pB) such 
that ipi G C and 

p p p 

^ (Pi(Xj) - (po^^X^j > - y^Wi(Xi). 
i=l i=l i=l 

It will also be convenient to consider 

(Q) inf | J (if)) , ip = (^l) •••>'0p) : eacn "0j is sucn that ipi — to, is convex j 
with 

r. P P 

v i=l i=l 77 

Note that by construction J(ip) > for every admissible ^ and J^) = if 
and only if 70 = 7^. 

Lemma 6.1. T7ie following holds 

r p 

max(P*) = inf(P) = inf(Q) - / (2^70(2;). 

Proof. Let us write (P) in the form 

inf F(A<p) + GUp) 

<p=(<pi,...,ip p ,ipo)eC(B)PxC(pB) 

where A : C(B) P x C(pB) — >■ C(B P ) is the linear continuous map defined by 

p p 
A(^(x) := ^^(xi) - y f ^ Xj j , Vx = (xi, G B p , 

i=l i=l 

and F, G are the convex lower semicontinuous (for the uniform norm) func- 
tional defined by 

F{6) = { J*^° if ,WeC(B>) 
[ +00 otherwise 
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It is easy to see that the assumptions of Fenchel-Rockafellar's duality theorem 
(see for instance [IS]) are satisfied and thus 



inf(P) = max -F*( 7o - 7) - G*(A*( 7 - To))- 

■y£M(BP) 

The adjoint of A, A* is easily computed as : M(B P ) ->■ M(B) P x M(pB) 
(where A4. denotes the space of Radon measures): 

A* 7 = (7 1 , 7 P , -n s7 ), V 7 G M(B P ). 

Direct computations give 

F*(7- 7o ) = { ~ 1*1*1 X 7>0 

+00 otherwise 

and 

Oif jeK{ lo ) 



G*(A*(7- 7o )) 



-00 otherwise. 



Therefore (V*) is the dual of (V) in the usual sense of convex programming 
and max(7 7 *) = inf('P). To prove that 

r p 

inf(P) = inf(Q) - / S^wiix^d^x), 
Jbp i=1 

take if G E and ipi := Wi + ipi for i = l, ..,p, the constraint then reads as 

P P 
y2^i( x i) > Pol^Xi), Vx G B p . 



i=l i=l 

Now in (V), one needs to choose <£> as large as possible without violating 
this constraint. Thus the best (po given (<pi, ...,<p p ) is 

which proves the desired identity. 

□ 

Lemma 6.2. Let ipi be such that ipi — Wi G C for every i and g = (g±, g p ) G 
C p . Then 



lim hjy, + Sg) - J(n = I 9i{*i)d{li ~ 7$) 

V 

^^(xi)d(7o -7^)(x). 



BP i= i 
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Proof. For 6 > 0, one first gets that 

hjy, + Sg ) - J(n = £ / ft(x0d(7S)- 

+ %)(x) - □i^(i))dm (a). 
pi? o v ' 

Note that the integrand in the second term is bounded since g is. Now 
fix some (xi,...,x p ) E B p , and set x = Y^,=\ x ii Vi := ^v>( x ) an d ^ 5 : = 
T^, + 5 g (x). Since Y^l^iVi = Y7i=iVi = x ' ^ comes as a direct consequence of 
the definition of infimal convolutions that: 

1 P 

-(□.(^ + S 9i ){x) - D^x)) < J^giiVi) (6.2) 

j=i 

and 

1 p 

+ - n^ixY) > Y,9i{vt)- (6.3) 

i=i 

Using the compactness of B and the strict convexity of ipi, it is easy to check 
that yf — > yi as 5 — > + . Therefore, from (16. 2 p and (16. 3ft one has 

1 P 

£m -(Di^i + <^)(x) - D t ^(x)) = $>(7$(x)) 



i=l 



and this holds for every x E pi?. It then follows from Lebesgue's dominated 
convergence theorem that 

1 P /" p r 

lim -[Jty + ^) - J{n = V / g^dinl) - V / ft (7j(x))dm (x) 

= 9i(xi)d(lo - = / $^(2^(70 -7^)0)- 

<=i ■ / * p i=i 

□ 

It follows from Lemma 16.21 that, if ip solves (Q), then 7^ dominates 70. 
Hence, if one knew that (Q) possesses solutions, the existence of an w-strictly 
comonotone allocation dominating 70 would directly follow. Unfortunately, 
it is not necessarily the case that the infimum in (Q) is attained - or at least 
we haven't been able to prove without additional conditions. The difficulty 
here comes from the fact that minimizing sequences need not be bounded (see 
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paragraph 14.41) . It may be the case that additional regularity assumptions 
on 7q would guarantee existence. In the present paper no such assumption is 
made, and a different path is chosen to overcome the difficulty by an appeal 
to Ekeland's variational principle. 

Lemma 6.3. Letting e > 0, there exists ip £ admissible for (Q) such that 

1. J(i> £ ) <mi(Q)+e 
2. 

r p 

\ipi{xi)d(^ e - 70 ) < 



limsup / y. 

£->0+ J BP i=1 



for every (ip 1 , tp p ) G C p 

3. 

f P 

liminf / V^^)^)^ - 7o) > 
£^o+ J BP 

for iff = ipi^ £ — Wi (these are convex functions by definition). 
Proof. For e > 0, let f e be admissible for (Q) and such that 

J(/ £ )<inf(Q)+ £ . 
Let then k £ > be such that 

lim i ek e [l+ \\f e \\] = (for instance k £ = - ). (6.4) 

e ^o+ e 7 C 1 + 11/611) 

It follows from Ekeland's variational principle (see [I7J and (5]) that for every 
e > 0, there is some ip £ admissible for (Q) such that 

He ~ fe\\ < p We) < J(fe) < mf(Q) + 8, (6.5) 

where \\h\\ stands for the sum of the uniform norms of the hi functions, and 
J{ip) > J{ipe) - k e £\\ip - ipe\\, W> = (ipi, ...,i/}p) : ipi - Wi e C, Wi. (6.6) 

Taking ip = ijj £ + 5ip with S > and tp G C p in ( 16. 6p . dividing by 5 and letting 
<5 — > + , one thus gets by the virtues of Lemma [6.21 



f P 

/ Y] <Pi(Xi)d(y - 7^) > -k t 

J BP i=1 



e\\ip\\. (6.7) 
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Using (16 .4p and letting e — > + yields 

f P 

limsup / y2<Pi(xi)d(^ e - 70) < (6.8) 

e^0+ J BP ~~J 

for every (tpi, ...,tp p ) G C p . To prove the last assertion of the lemma, write 
ip £ = (p e +w with Lp e G C p . Then for 5 G (0, 1) one has ip e —5ip £ = (l—5)(p e +w, 
and then (16.61) can be applied to ip e — 5(p e , leading to 

l -[J(i, £ -5^)- J(V0] > -Ml¥>1l> 

which, letting 5 — > + and using the same argument as in lemma 16.21 leads 
in turn to 

f P 

By ([621) and (IB31) . it follows that 

fc B e||^|| < Mlkll + ll^-MI + IIMI) < MIHI+e+MlAll -+ as £ + , 



which enables us to conclude that 

p 



liminf / V (pl(xi)d(^ s - 70) > 0. 
e ^o+ J BP 



(6.9) 



□ 



Lemma 6.4. Let ip e be as in lemma ROI and set 7 £ := 7^ then up to some 
subsequence, 7 e weakly star converges to some 7 (w-comonotone by construc- 
tion) such that 7 G .Ms (mo) and 7 dominates 70. Moreover 7 solves (V*). 

Proof. By the Banach-Alaoglu-Bourbaki theorem, one may indeed assume 
that 7 £ weakly star converges to some 7. Obviously, 7 is i/7-comonotone and 
IIs7 = Ils7o = m , hence 7 G .Ms (mo). The fact that 7 dominates 70 is 
obtained by letting e — > + in (16.81) . It remains to prove that 7 solves (V*). 
Defining ip e : = ip £ — w as in Lemma 16.31 one has 



/ V^(xi)d(7o-7e) + / rjd(j e -To) inf(Q) as £ ->• 0~ 
By (16. 9p . passing to the limit thus yields 

inf(Q) < / ^(T-7o), 

J BP 
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which, combined with Lemma 16.14 gives 



r]dq > inf(Q) + / 77^70 = max (P*). 
bp J bp 



□ 



Lemma 6.5. Let 7 be as in lemma (TJ. Then: 

1. if '70 solves (V*) then 70 is w-comonotone, 

2. 7 strictly dominates 70 unless 70 is itself w-comonotone. 

Proof. If 70 solves (V*), it follows from Lemma I6TT1 that inf(Q) = 0. For any 
minimizing sequence ip £ (not necessarily the one constructed in Lemma [6. 8 1) 
of (<2), the following holds: 

C v p 
= lim J(ip e ) = lim / ( V^ je (x;) - ) )^7o(^) 

8=1 8=1 

= s 5? + / BI (t^)-E^(n,(E-<))) 

J B ' 8=1 8=1 8=1 



P 

d7o(x). 



By a density argument, one may consider ip £ a minimizing sequence such 
that each ip £ belongs to C 1 (5). Fix (xi, ...,x p ) and set x := Y^=i x i- Then 
y e := T^ e (x) can be characterized by the fact that there is a vector p £ R d 
and a vector (Aj) of nonnegative weights such that 

p 

V^, £ (yf) = p - A^ e , Aj = if yi $ dB, J^vt = x - ( 6 - 10 ) 

8=1 

On the other hand, since Wi is strictly convex and tp ii£ — W{ £ C, it follows 
that for any a and b in B 2 , 

if> i>E (b) - if> i>£ (a) > V^, e (a) • (b - a) + 9i{\b - a\) (6.11) 

where function Bi is defined by 

Biit) := mf{wi(b) - uii(a) - Vwi(a) • (b - a), (a, b) £ B 2 , \a - b\ > t} 

for any t £ [0, diam(B)]. Function Bi (the modulus of strict convexity of 
Wi) is a nondecreasing function such that Bi(0) = and 9i(t) > for t > 0. 
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Combining (jfi.lUp and (16. lip yields 

V P V V 

i=l i=l i=l i=l 

VP V 

i=l i=l i=l 

>X>(|*,-^|). 
i=l 

Hence the fact that J{ip £ ) -> as e -> + implies 

Yo{x) = 



J BP l=l j 



so that 

TfsC^j x j) — ^ — ^ as e — )■ + for 7 -a.e. x. 

j 

Therefore, by Lebesgue's dominated convergence theorem, 

f(x)dy (x)= lim / f(T^(S^Xj))d^ (x) = lim / f(T i)e (x))dm (x) 

lim / fd-f A 



Bp 



holds for all / G C(B P ). Hence, 7^ e weakly star converges to 70 which proves 
that 70 is w-comonotone and Point 1 is proven. We now prove Point 2. If 
70 is not w-comonotone then by Point 1, it does not solve (V*) and thus 
J rjdipf — 70) > so that 

r P r P 

/ z2 w i{ x i) dn f < / ^ J Wi(Xi)d'Jo 

Jbp i=1 Jbp i=1 

hence 7 strictly dominates 70. This completes the proof. 

□ 



6.4 Proof of Theorem 13.1 

Let 70 := £(X) = C(X\, X p ) and consider again the minimization problem 
(V*). It follows from Lemma [6.41 that there exists a w;-comonotone solution 
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7 to (V*). By construction, 7 dominates 70, and by the definition of w- 
comontonicity, there is a sequence 7 n of u>-strictly comonotone allocations 
that weakly star converges to 7. Lemma 14.11 allows to write 7™ = C(Y n ) 
for some Y n G A(X) which is obviously comonotone in the univariate sense 
according to Definition 12.41 Since Y n is bounded in L°° and each Y n is a 1- 
Lipschitz function of X, one may assume that Y n converges uniformly (up to 
a subsequence) to some Y 6 A(X). It is obvious that Y is also comonotone 
and that 7 = £(Y), hence Y dominates X. The strict dominance assertion 
follows from Lemma [6.51 
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